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Two recent reports confirm that the newly synthesized Sr2Co04 (formal oxidation state Co*"*") 
shows a high Curie temperature (^ 250 K), but they report different moments of 1.8 and 1 
per Co. Using both commonly used functionals in the correlated band approach (LDA+U) as well as 
the local density approximation (LDA), the combined effects of correlation and hybridization with 
O 2p states are calculated and analyzed. Sr2Co02 is already ferromagnetic within LDA (Af=1.95 
/is). Increasing U from zero, the two LDA+U schemes affect the moment oppositely out to a 
critical value (7c=2.5 eV, at which point they transform discontinuously from different states to the 
same large U state. Fixing U at f/c, fixed spin moment calculations show similar behavior out to a 
minimum at I^lb (a half metallic state), beyond which the fuUy-localized-limit scheme jumps to a 
state with energy minimum very near 2/is (very close to the LDA moment). Although the energy 
minima occur very near integer values of the moment/Co 2/i_B), the strong 3d — 2p mixing 

and resulting 3ci orbital occupations seem to preclude any meaningful 5 = 5 or 5 = 1 assignment 
to the Co ion. 

PACS numbers: 71.20.Bo, 71.27.+a, 71.30.+h, 75.30.Gw 



I. INTRODUCTION 

First row (3(i) transition metal ions in high oxida- 
tion states have been of interest for some time, due 
to their competing and delicate spin states. There 
are still several cases where behavior is not under- 
stood, indeed sometimes (due to sample questions 
arising from difficulty in synthesis) the data is not 
unambiguous. The best known, and most thor- 
oughly studied, example is that of the quasi-two- 
dimensional (2D) cuprates, which when hole-doped 
become high temperature superconductors. In their 
undoped state they are Cu^"'"-based antiferromag- 
netic insulators. Hole doping drives the Cu oxidation 
state toward the unstable (practically non-existent) 
Cv?+ state. At a doping level of 0.15-0.20 holes/Cu 
optimal superconductivity is reached. They can be 
doped beyond that level, when they become con- 
ventional metallic Fermi liquids. Upon hole doping 
the 'oxidation state' designation must be interpreted 
with care, since it is clear that the holes go onto the 
oxygen ions to a large degree. 

Another high-oxidation-state ion is 00"*+. This 
ion has recently gained wide exposure due to 
the unusual properties of Na2;Co02, including the 
discovery Q that this system becomes superconduct- 
ing (4.5 K) around a:=0.3 when it is hydrated. 
This value of x suggests that the Co ion is 70% 
of the way to being Co'*+, or alternatively, 70% 
of the ions are 00*+ while 30% are 00^+. This 
system is however metallic for all x except for 
Na ion/orbital/spin/charge orderingQ precisely at 
a;=0.5 (precisely what is responsible is not yet clear). 
There is magnetism and correlated electron behav- 



ior for X > 0.5 (the Co-^^ end) while for x < 0.5 the 
materials appear to be weakly correlated nonmag- 
netic metals. While there has been much expecta- 
tion that the (metastable) endpoint member C0O2 
(nominally Co^+jd^) is a Mott insulator, the evi- 
dence is that it remains a nonmagnetic metal. [S] Cal- 
culations indicateQl that, as in the cuprates, upon 
hole doping from the Co^+ end, much of the charge 
difference occurs on the O ion. 

The La2-a;Sra;Co04 system (LSCoO) with two- 
dimensional layered K2NiF4 structure has become 
interesting because of its magnetic and electrical 
properties. In this system the formal oxidation 
state Co^^"'"^-'+ near the La end is not unusually 
high, so its behavior might be expected to be read- 
ily understandable. With increasing Sr concentra- 
tion X, LSCoO shows a structural transition from 
orthorhombic to tetragonal (i.e., a lattice constant 
6— >a) around x=0.5, with enhanced two-dimensional 
electronic properties. 5] The structural transi- 
tion may be connected to an antiferromagnetic- 
ferromagnetic transition (the maximum Curie tem- 
perature rc=220 K at a;=0.9),@ accompanying a 
magnetic change that has been interpreted in terms 
of a Co^+ spin-state transition from high spin to 
intermediate spin configurations around j:=O.7.0 
There is another suggestion that at x=l the sys- 
tem is a high-spin, low-spin charge-ordered state. Q 
However, there is no agreement amongst the mea- 
surements on the magnetic behavior |^ 0, ITU 
nor metallic behavior|6l Il2|. The end member x=0 
is an antiferromagnetic (AFM) insulator with Neel 
temperature T^—^Tb K. 13] 

Although earlier studies were confined to the 



range below 2:= 1.4, recently the end member 
Sr2Co04, formally Co'^"*", was synthesized by Mat- 
suno et al.^^ [I3 and by Wang et al.^^ Iit\ . 
and characterized as FM with high Curie tempera- 
ture Tc~250 K. However, their differing synthetic 
methods have led to different properties. Mat- 
suno et al. synthesized a single-crystalline thin 
film using pulsed-laser deposition, while Wang et 
al. produced polycrystalline samples under high 
pressure, high temperature conditions. The former 
shows metallic T-dependent resistivity below Tc, al- 
though it has definitely higher resistivity of order 
of 10^^ — 10""^ cm at low T than in a typi- 
cal metal. The pressure-synthesized samples show 
nearly temperature independent resistivity (perhaps 
due to polycrystallinity), but with similar magni- 
tude. The more confusing difference is in the ob- 
served saturation magnetic moment, 1.8 /zb and 1 
Hb/Co respectively. Nevertheless, the sample with 
the smaller ordered moment has been observed|l6j 
effective (Curie- Weiss) moment Pe//= 3.72 //s, char- 
acteristic of a much higher S = 3/2 spin configura- 
tion. Such a moment would suggest a much higher 
ordered moment, < Sx >~ 3^ib- 

The related perovskite system Lai-a-Sr^rCoOa 
(formally Co(3+^)+) has been studied for some 
time.,18„ Jl, JHH |2^ 23, magnetic prop- 

erties are also altered by Sr concentration x: non- 
magnetic for X < 0.05, spin glass for 0.05 < x < 0.2, 
and ferromagnetic above a;=0.2. For the end mem- 
ber, SrCoOs ("Co'*'*'") shows metallic conductivity 
and has magnetic moments 1.25 and 0.1 hb for Co 
and O respectively, 13 whereas the ground state of 
LaCoOs is a nonmagnetic band insulator. Its exci- 
tations have been interpreted in terms of a locally 
orbitally-ordered excited state. p^ l25l | 

In this paper, we look in some detail at the elec- 
tronic and magnetic structure of Sr2Co04 both as 
uncorrelated, using the local density approximation 
(LDA), and viewing effects of correlation as de- 
scribed by the LDA-I- Hubbard U (LDA+U) method. 
Both commonly used LDA-I-U functionals are em- 
ployed, and their results are compared and con- 
trasted. A number of unusually rich magnetic phe- 
nomena arise (including a half metallic phase), re- 
flecting the strong hybridization with O 2p states 
that complicates the accommodation of correlation 
effects on the Co ion. 



II. STRUCTURE AND CALCULATION 

Sr2Co04 has the bet K2NiF4-type structure, 
space group lA/mmm (No. 139), pictured in Fig. 
□ The Co - planar O (PO) bond length (1.878 A) 
is shorter by about 6% than that of Co - apical O 




FIG. 1; (Color online) Crystal structure of Sr2Co04, a 
body-centered-tetragonal type. The planar 0-Co bond 
length is 1.878 A, about 6% shorter than the apical 
O-Co bond length. Lattice constants used here are 
a=3.755 A and c=12.6 A. 

(AO) (1.988 A). The distorted CoOe octahedron, 
elongated along c-axis, leads to crystal field split- 
ting of t2g Eg{dxz,dyz) + B2g{dxy) aud eg -> 
Ag{d^z^_r2) + Big(d^2_j^2) states. We have used the 
lattice constants a=3.755 A, c=12.6 A, and apical 
O (0.1578) and Sr (0.3544) internal parameters op- 
timized by Matsuno et aZ.,|Q| which are consistent 
with the experimentally measured values by Wang 
et a/, [ill 

Our calculations were carried out within the lo- 
cal density approximation (LDA) and LDA-I-U ap- 
proaches with the full-potential nonorthogonal local- 
orbital (FPLO) method. "^G*! Both commonly used 
schemes of the LDA-f U method were employed so 
comparisons of the predictions could be made. Both 
forms have the same Hubbard-like density-density 
interaction 

7) ^ ^ \U^uin' ^(T,cr' <-^mm']^mcr^m'cr'(il) 

where {rima} is a site occupation set and J is in- 
tra exchange integral. The form displayed here 
is schematic in the sense that it does not display 
all of the indices involved in the full coordinate- 
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system-independent form that is implemented in the 
code. The two LDA+U approaches differ only in the 
method of treatment of the "double counting" term, 
intended to subtract out the shell-averaged interac- 
tion that has already been included in LDA. One 
choice is the so-called "around mean field" (AMF) 
scheme, which is expected to be more suitable when 
the on-site Coulomb repulsion U is not so strong. j23| 
The other choice, called the "fully localized limit" 
(FLL) (also called "atomic limit"), is more appro- 
priate for large U systems, psj The methods to treat 
the double counting problem in the both schemes 
are given by 
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where n is the shell-averaged occupation, while n„ is 
its spin-decomposed analog. These double-counting 
terms can be written in other forms that empha- 
size other aspects of the interaction. However, 
this form is illustrative because it emphasizes that 
the difference lies in the magnitude of the (self- 
consistent) atomic moment. Clearly it is only the 
spherically averaged values of U and J that enter 
the double-counting terms. The difference between 
the two forms is that the double-counting term, and 
the resulting potential, includes a spin dependence 
in the FLL form. 

The fully relativistic scheme of FPLO'SO^ was 
used when spin-orbit coupling (SOC) in LDA 
was necessary. The choice of basis orbitals 
were Sr (4s4p)5s5p4(i, Co (3s3p)4s4p3(i, and O 
2s2p(3s3p3(i). (The orbitals in parentheses indicate 
semi-core or polarization orbitals.) The Brillouin 
zone was sampled by a regular mesh containing up 
to 641 irreducible k points for LDA and LDA-fU, 
and 1639 irreducible k points for SOC and fixed spin 
moment (FSM)|2l| calculations that require more 
carefully treatment near the Fermi level {Ep)- 
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FIG. 2: (Color online) FM LDA majority (upper) and 
minority (lower) band structures along symmetry direc- 
tions. The thickened (and colored) lines emphasize Co 
dzz'^-r^ (green or light) and d^i_y2 (red or black) char- 
acters, which form bonding and antibonding bands with 
apical and in-plane O Pa states respectively. The sym- 
metry points for the body-centered-tetragonal structure 
follow the Bradley and Cracknell notation as given in 
Fig. |H] The dashed horizontal line denotes the Fermi 
energy. 



III. UNCORRELATED TREATMENT 

A. Electronic structure 

The FM LDA band structures, exhibiting total 
magnetic moment A/=1.95 /is that is only acciden- 
tally near an integer value, are shown in Fig. [2| 
The five Co 3d bands and three 2p bands of each O 
constitute an entanglement of 17 hybridized bands 
for each spin, with total p — d bandwidths of 10 eV 
and 11 eV for the majority and minority bands re- 
spectively. The strong hybridization throughout the 



bands is evident in the corresponding total DOS and 
accompanying atom-projected DOS, shown in Fig. 
The crystal field splitting of 2 eV, identified from 
the density of states, between and Cg manifolds 
is only a little higher than the exchange splitting 1.3 
cV. As noted earlier, these t2g and eg designations 
are broken down by the 4/m symmetry of the Co 
site. The electronic structure shows clear quasi-two- 
dimensionality, consistent with resistivity measure- 
ment by Matsuno et al. 

In Fig. [51 the pda and pda* bands (for each of 
d^2_y2 and d3^2_^2 orbitals) are highlighted. The 
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FIG. 3: (Color online) Total and atom-projected densi- 
ties of states for LDA FM calculation. Near Ef, there 
is a van Hove singularity in the minority channel. While 
AO is almost fully occupied, the minority of PO is par- 
tially occupied, resulting in large magnetic moment for 
PO (for details, see text). Additionally, the band width 
of AO is by 40% less than that of PO. The DOS at Ef 
N(0) is 2.90 states/eV per both spins. 



more dispersive bands (darker, or red) arise from 
bonding and antibonding interactions between Co 
dx^-y^ and PO pa states. The bonding and anti- 
bonding bands have a separation of 10 eV at the X 
point, this total width arising from a combination of 
the differences between the t2g and 2p site energies, 
and the hybridization between them. The dispersion 
of the antibonding (upper) band can be described 
roughly by effective hopping amplitude t = 0.40 eV 
(majority) and t — 0.53 eV (minority). The behav- 
ior of the (1^2 _y2 band is similar to that observed 
in the high superconductor La2Cu04, in which 
this antibonding band plays a central role,[33, 113 
but it is mostly unoccupied here. However, with 
the lower number of 3d electrons the cobaltate elec- 
tronic structure is influenced more strongly by the 
t2g manifold, which forms relatively flat bands, the 
top of which hovers around the Fermi energy. In the 
minority bands a van Hove singularity lies at Ep at 
the Z point, perhaps contributing to the positioning 
of the Fermi level and hence to the total moment. 

The less dispersive highlighted bands (green, or 
lighter), lying in the - 2 eV range, arise from the 
antibonding interaction of Co (i3^2_^2 and AO 
states. The bonding band of this pair consists of a 
remarkably flat band at -4 eV for each spin direction. 
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FIG. 4: (Color online) Orbital-projected densities of 
states for Co 3d (upper) and O p states (lower) in LDA 
FM calculation. The d^cz and dyz states are degenerate. 
The crystal field splitting of 2 eV between t2g and eg 
manifolds is only a little higher than the exchange split- 
ting 1.3 eV. The PO px and AO pz are the cr-orbital. 



B. Magnetic Tendency 

The FM state is favored energetically over the 
nonmagnetic state by 0.37 eV/Co, similar to the 
value given by Matsuno et al. The spin magnetic 
moments are 1.95 total, 1.52 from Co and 0.22 from 
each PO (in units of hb)- The large magnetic mo- 
ment for PO, due to strong hybridization with Co 3d 
bands, has been observed previously in LiTCu09[3^ 
and a few other cuprates. The MuUikan decomposi- 
tion of O charges as well as the band filling of related 
bands indicates that AO is consistent with its for- 
mal designation 0^~, whereas PO contributes con- 
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TABLE L Co 3d orbital (MuUikan) occupancy in LDA, 
where Af=1.95 fiB- The difference of occupancies be- 
tween both spin channels is directly related with contri- 
bution of each orbital to spin magnetic moment, which 
can be seen to be spread over all five 3d orbitals. 
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FIG. 6: (Color online) LDA FM Fermi surfaces, for 
M= 1.95 fj,B- Surfaces (a) and (b) are from the ma- 
jority states, surfaces (c)-(f) from the minority bands. 
While (a) and (c) contain holes, the others enclose elec- 
trons. Each surface has mainly (a) d^2_y2 (0.5 holes), 
(b) d^^2_j.2 (0.01 electrons), (c) dxz (0.2 holes), (d) d^y 
(0.4 electrons), (e) dy^ (0.4 electrons), and (f) d^2_y2 
(0.02 electrons) characters. The number in parentheses 
says carrier number containing each Fermi surface per 
Co. 



FIG. 5: (Color online) Fixed spin moment calculations 
in LDA and both LDA+U schemes at C/c=2.5 eV. The 
arrows pinpoint (meta)stable states. The dashed line for 
LDA FSM indicates fitting line with E{M) - £(0) = 
eo - ctM'^ -I- /JM", where £o=5 meV, q:=133 meV/^i|, 
and /9=16 meV/zig. 



siderably to the conduction bands and cannot be 
considered fully ionic. According to the Co 3(i or- 
bital occupancies given in Tabled however, every 3c? 
orbital contributes to the magnetic moment due to 
itineracy. The variation from average contribution 
occurs in two states: dxy has 40% larger, and d'^^2_^2 
orbital 40% smaller, contributions than the average 
contribution to the moment. The strong itinerant 
character (all Sd orbitals are neither fully occupied 
nor fully unoccupied) explains why no Jahn- Teller 
distortion is observed in Sr2Co04. 

The fixed spin moment method[3H is applied in 
a following subsection to probe the magnetic be- 
havior. In the LDA fixed spin moment calculations 
the low moment region is given by E{M) — E{0) w 
-aM^+PM^ with constants a=133 meV//x|, ^=16 



meV//z^. The Stoner-enhanced susceptibility is 
given by 

X = Xo/[l - iV(0)/] (3) 

where the bare susceptibility xo = 2/x^A^(0) and 
A'^(O) is the single-spin density of states at the Fermi 
level. Also in the low-M limit one obtains formally 
a = (l/2)x~^, giving the Stoner interaction / = 
1.17(±0.05) eV. With this value /iV(0)« 1.7, giving 
a very strong Stoner instability of the nonmagnetic 
phase that is numerically similar to that of nonmag- 
netic Fe. 



C. Fermiology 

The LDA FM Fermi surfaces (FS) pictured in Fig. 
|S1 consist of two sheets from the majority states and 
four minority sheets, with very simple geometry and 
strong two-dimensionality. Except for one sheet that 
has an ellipsoidal shape (majority d^^^-r^ charac- 
ter), the FSs have the shape of rectangular cylin- 
ders with rounded corners. The X-centered sheets 
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FIG. 7: (Color online) LDA FM band structures from 
spin-orbit coupling calculation with the quantization 
axis along (100) (solid) and (001) (dashed) directions. 
Note the (001) direction is the easy axis. The dashed 
horizontal line denotes the Fermi energy. 



contain holes, whereas the F-centered surfaces con- 
tain electrons. The FS arising from majority dj.2_y2 
(Fig. Ela)) has a clear nesting feature, as does the 
minority hole sheet in Fig. Eld). The spanning vec- 
tors may lead to spin-density- wave or/ and charge- 
density-wave instabilities. Intra-surface scattering 
may also show some nesting features. For the most 
part Fermi velocities vp are in the range of a few 10'' 
cm/sec as usual in a metal, but order of magnitude 
lower velocities occur along the F3 to Z line at Ei?. 



D. Magnetocrystalline anisotropy 

We included spin-orbit coupling by using the fully 
relativistic FFLO method 30] to allow us to study 
the strong magnetic anisotropy reported^^ by Mat- 
suno et al. The total energies along (100) and (001) 
directions were calculated to determine the magne- 
tocrystalline anisotropy energy (MAE). The energy 
difference _E^ioo) ~'E'{ooi>=0-57 meV/Co is consistent 
with experimentally observed (001) easy c-axis. The 
relativistic calculations give an orbital magnetic mo- 
ment of 0.08 /is for Co, and negligible values for oxy- 
gen ions. The spin and orbital magnetic moments 
depend on the spin direction, differing by about 5%. 

The experimental MAE can be estimated from 
field dependent magnetization along the two 
directions. Tl'l The data show half-saturation at 7 T 
along the (100) direction. Considering the saturated 
magnetization to be Ms~1.8 /is/Co, the MAE from 
Mg-B is 1.5 meV/Co, a value three times larger than 
our calculated result. Such underestimation in LDA 
is common for Co compounds. The reason is still 
unsettled, but candidate explanations are the poor 
treatment of orbital Dolarization^35.J or for Hubbard- 
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FIG. 8: (Color online) Effect of the on-site Coulomb 
repulsion U on total and Co local magnetic moments 
in both LDA-I-U schemes. At Uc=2.5 eV, a metal to 
half metal transition occurs. The first-order transition is 
obtained nearly at the same U in the both schemes. (In 
fact, the transition occurs a httle higher Uc in FLL, but 
the difference is only less than a few tenth eV.) 



like correlation in the standard LDA. 

Rotating the spin direction induces a change of 
band structure, as shown Fig. [7| While the band 
structure for the moment along (100) direction is es- 
sentially that of LDA, there are additional splittings 
at several points for the moment oriented along the 
easy axis. The most visible splitting occurs in Co 
t2g manifold, where the spin-orbit splitting spans 
Ei? at the Z point. Applying SOC with spin di- 
rected along the easy axis, the bands split into pri- 
marily (ij=5/2,mj=-3/2 character in the upper band, 
and dj^3/2,mj=i/2 in the lower band. The splittings 
are 60 meV and 90 meV at the Z and F point re- 
spectively. 



IV. INCLUSION OF CORRELATION 
EFFECT 

A. Metal to Half-metal Transition 

Since the appropriate value of U in this and other 
cobaltate systems is unclear, we have studied the 
ground state as U is increased. The moment initially 
increases slightly in FLL, while it decreases slowly to 
1.5 HB at C/=2.5 eV in AMF, as shown in Fig. |H1 
From their separate states just below and at the crit- 
ical value Uc=2.5 eV, beyond this critical value the 
moment drops sharply to 1 /xs in both schemes, and 
then both schemes produce very similar results in 
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FIG. 9: (Color online) LDA+U FM majority (upper) 
and minority (lower) band structures at (7c=2.5 eV in 
the half metallic state with M=l ^b- The majority state 
shows a gap of 0.25 eV. The thickened (and colored) lines 
emphasize Co d2,^-2_^2 (green or light) and d^2_y2 (red 
or black) characters. The dashed horizontal line denotes 
the Fermi energy. 



the entire region above Uc- Note particularly that 
this transition starts from distinct states, but oc- 
curs at the same value Uc to the same final state. 
Both high-moment and low-moment states can be 
stabilized in the calculations at Uc- Unlike studies 
in the Naa;Co02 system, no discernible hysteretic 
region could be found at this first-order transition. 
This magnetic collapse accompanies a metal to half 
metal transition, presumably because there is a par- 
ticular stability of this half metallic (HM) FM state 
since both LDA-I-U schemes transition to it. The 
HM state, with a magnetic moment 1 /^s, has been 
observed also by Wang et al. As U increases, the 
Co local magnetic moment decreases to vanishingly 



TABLE II: Co 3d orbital occupancy in half metallic 
state, having M=l fiB, at U—2.5 eV. Compared with 
LDA result shown in Table |I] there are two remarkable 
changes in the eg state; nearly vanishing contribution 
of d^^2_j.2 and negative contribution of d^2_y2 (for de- 
tails, see text). It seems to be close to LS state, but 
the d^2_y2 minority has considerable occupancy due to 
strong hybridization (itinerary) which makes impossible 
to be called strictly as LS state. The total occupation is 
6.79. 
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small value by U ^6-7 eV. The state remains a HM 
FM, the moment has been pushed onto the PO ions. 
Since the moment on Co has vanished, it is not sur- 
prising that the two LDA+U schemes give the same 
result. 

The microscopic mechanism behind the magnetic 
collapse induced by the on-site Coulomb repulsion 
U can be unraveled from a study of the charge de- 
compositions in Tables and ^ and comparison 
of band structures. We can compare and contrast 
the two viewpoints. First we point out that, while 
the MuUikan charges given in these tables are some- 
what basis set dependent so their specific magni- 
tude should not be given undue significance, differ- 
ences - whether between orbitals or between spin- 
projections - are more physical. Then, while these 
orbital occupations provide one characterization, we 
have also provided in Figs. [21and|51by the fatbands 
technique, the bands that one identifies with dx2-y2 
and d^^2_^2 character. Using both viewpoints pro- 
vides a more robust interpretation of behavior than 
either separately. 

From the tables one can determine that increasing 
U does not change the total 3c? occupation, but intro- 
duces rearrangements in the 3d charge and moment. 
From the bands, the main difference (between FM 
LDA and HM FM LDA-hU) is that the dispersive 
majority d^2_y2 band has become fully MTTOCCupied, 
leaving the small gap that results in half metallic- 
ity. Returning to the charges and moments, what 
stands out is that the moment on the d^^^y'^ or- 
bital has flipped, from -f0.35^B to -0.43/iB. The 
c?3^2_^2 moment has undergone a smaller change in 
the same manner: to -0.09/l<_b. Thus the net 

moment on the Co ion, which is ~ 0.5^b, is the re- 
sult of ~ +1^J,B in the t2g orbitals (primarily dxy) 
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and ~ — 0.5/iB in the eg orbitals. This type of can- 
cellation, in more striking form, has been seen previ- 
ously in LDA-I-U results for LaNi09.[37l| which has 
an unusually low formal oxidation state for a nicke- 
late. 



B. Fixed Spin Moment Calculation at Uc 

The results of the previous section show that two 
magnetic states coexist at Uc in LDA-I-U, whether 
one uses the AMF or FLL functional. Such a coexis- 
tence has been already observed in LDA-I-U calcula- 
tions for the sodium cobaltates, where the change 
at Uc corresponds to a charge disproportionation 
transition. 38] Here the change is simply in the state 
of the (single) Co ion; there is no experimental indi- 
cation of disproportionation here. Here we analyze 
the FSM results at Uc for LDA, AMF, and FLL. 

The energy versus total magnetic moment behav- 
ior, displayed in Fig. [SJ shows very interesting dif- 
ferences as well as similarities. Perhaps most inter- 
esting is that positions of local minima occur at (or 
near) integer values M=0, 1, and 2 (in units of /is), 
as if there might be underlying states with Sz — 0, 
i, or 1 pervading the behavior. The results of previ- 
ous sections however established that there is strong 
d — p hybridization which renders integral moments 
no more favored than other values. The LDA curve 
shows simple behavior: a Stoner instability at M=0 
and a single minimum near (but not precisely at) 
M=2. The LDA-I-U results show more interesting 
variation. 

For both AMF and FLL the paramagnetic state 
is metastable, with energy rising up to a moment 
M X 0.5, and then decreasing very similarly to the 
minimum at M=l. Beyond that point, AMF levels 
off to a broad flat region M=1.3-1.8 beyond which 
it increases with M . FLL however switches over to 
a separate phase, with its minimum at M=2. The 
E(M) curve for FLL seems to consist of two sepa- 
rate parabolas (different phases) with minima very 
near M=l and M=2. The M=l result is integral be- 
cause a half metallic phase is encountered (see Fig. 
inj and not due to an 5* = 1/2 configuration of the 
Co ion. Both majority dj;2_yi and d-^z^_r-2 become 
completely unoccupied, leaving a small gap to the 
t2g bands. 

We must point out that for both AMF and FLL 
schemes, the minimum at M—Q fiB does not corre- 
spond to a nonmagnetic state. There is a moment 
on Co with magnitude 0.14 fiB that is canceled by 
magnetization on the PO ion. The driving force for 
favoring this low, canceling moment phase over a 
nonmagnetic state is not clear. Possibly an anti- 
ferromagnetic result, requiring a doubled unit cell, 
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FIG. 10: (Color online) Effect of strength of U on the 
orbital-projected Co magnetic moment Aim, defined by 
difference between majority and minority occupancies, 
of Co d^2_y2 and dxy states in AMF scheme. (Results 
for the FLL scheme are similar.) In the large U limit an 
on-site "singlet" type cancellation of moments occurs in 
this m — ±2 channel (see Fig. |H]and text). 



would be a lower energy M=0 solution. 



C. Strong Interaction Regime 

In Sr2Co04 correlation effects should not be very 
strong it seems, since LDA already gives a FM state 
that seems consistent with one of the experimental 
reports. However, in a strong correlated regime, but 
not beyond realistic range of U, there is another in- 
teresting feature. Figure [TUI shows the JJ-dependent 
orbital- projected local magnetic moment Mm, de- 
fined by the difference between majority and minor- 
ity occupancies, of Co dx2_y2 and dxy states. Upon 
increasing U, the d^y minority state loses electrons, 
while the (i^2_j,2 minority state gains electrons. Be- 
yond U — Uc = 2.5 eV, dxy and d^^-y^ have large 
positive and negative local magnetic moments re- 
spectively, leading to an on-site "singlet" type can- 
cellation within this |m| — 2 channel. This type of 



cancellation, but within the Cg manifold, has been 
seen already in our previous results for LaNi09.[3^ 
In contrast with LaNi02, which shows the cancella- 
tion in the total magnetic moment and Ni^+^Ni^+ 
conversion, the moment of the system remains un- 
changed because of large magnetic moments on pla- 
nar oxygen ions. 
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V. SUMMARY 

Synthesis of the high formal oxidation state com- 
pound Sr2Co04, bulk materials by high-pressure 
high-temperature techniques and films by pulsed 
laser deposition (PLD), have led to a high Curie tem- 
perature metallic ferromagnet that introduces new 
transition metal oxide physics and may be useful in 
spin electronics devices. We have provided an in- 
depth study of the electronic and magnetic structure 
of this compound, looking specifically into the com- 
bined effects of correlation on the 3d orbitals and 
strong hybridization with O 2p states. 

Within LDA Sr2Co04 is metallic with a ferromag- 
netic moment near 2/iB, close to the saturation mag- 
netization reported for the PLD films. Application 
of the two commonly used LDA-I-U functionals re- 
veals several surprises. As U is increased from zero, 
the two functionals produce changes in the moment 
of opposite sign up to the critical value Uc=2.5 eV. 
Beyond half metallic ferromagnetic phase with mo- 
ment Within this phase, increasing the value of 
U has the effect of pushing the moment completely 
off the Co and onto the planar O ions by i7 ~ 6-7 
eV, while the total moment remains fixed at IfiB- 

With U fixed at J7c=2.5 eV, fixed spin moment 



calculations were carried out for both LDA-I-U func- 
tionals, and compared with the corresponding LDA 
result. Both LDA-I-U schemes behaved similarly out 
to a minimum at I^b (the half metallic state). Be- 
yond this the two functionals departed in their be- 
havior, with the FLL scheme jumping to a new state 
with minimum very near not half metallic and 
much the same as the LDA minimum (1.95/LtB). In 
the stronger-interacting regime U > 3 eV (which 
may not be appropriate for Sr2Co04), the corre- 
lation described by the LDA-I-U approach leads to 
oppositely directed moments on the dx-^-yi and dxy 
orbitals, reflecting the strong difference in hybridiza- 
tion of these two orbitals. 
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